Quantum electrodynamics corrections to energies, transition amplitudes and parity 
nonconservation in Rb, Cs, Ba + , Tl, Fr and Ra + . 
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We use previously developed radiative potential method to calculate quantum electrodynamic 
(QED) corrections to energy levels and electric dipole transition amplitudes for atoms which are 
used for the study of the parity non- conservation (PNC) in atoms. The QED shift in energies and 
dipole amplitudes leads to noticeable change in the PNC amplitudes. This study compliments the 
previously considered QED corrections to the weak matrix elements. We demonstrate that the QED 
corrections due to the change in energies and dipole matrix elements are comparable in value to 
those due to change in weak matrix elements. 

PACS numbers: 11.30.Er, 31.30.jg 



I. INTRODUCTION 

The study of parity nonconservation (PNC) in atoms 
provides a low-energy search for new physics beyond the 
standard model (see, e.g. [H, 0), which is a relatively 
inexpensive alternative to the high-energy searches per- 
formed in colliders. In fact, PNC in cesium is currently 
the most precise low-energy test of the electroweak the- 
ory due to the high accuracy of the measurements 0, 
and the calculations needed for their interpretation 
(see also @-||). 

The level of precision that has been obtained in atomic 
physics calculations and measurements has meant that 
strong-field QED corrections are observable and must be 
taken into account. It was shown in fact, that the inclu- 
sion of self-energy-type QED corrections to PNC calcula- 
tions in cesium restored the agreement with the standard 
model [I3-Q3 (see also [Pjl - fla ]). Just as these calcula- 
tions have proven to be important in the case of cesium, 
they will be necessary for the calculations of other atoms 
as the accuracy increases and as new experiments become 
available. For this reason, using the "radiative potential" 
method developed in Ref . , we present calculations of 
the QED corrections to the PNC amplitudes of several 
transitions in Rb, Cs, Ba + , Tl, Fr and Ra + . 

The case of rubidium is interesting because of its simple 
electron structure and small value of different corrections 
to the PNC amplitude @. The interpretation of the PNC 
measurements for Rb can be more reliable than for other 
atoms. We have shown in our previous work Q that 
the accuracy of the calculations for rubidium can surpass 
those for cesium, while the PNC amplitude is only several 
times smaller § (see also [Hj]). 

Francium is a particularly important application. The 
FrPNC collaboration has begun the construction of a 
laser cooling and trapping apparatus at the TRIUMF 
laboratories in Canada with the purpose of measuring 
atomic parity nonconservation in artificially produced 



francium [T3|. With a PNC amplitude expected to be 
around 15 times larger than that of cesium, and its rela- 
tively simple electronic configuration which leads to accu- 
rate calculations, francium is an ideal atom for precision 
measurements of PNC [18l42C| . When these measure- 
ments become available it will be very important to have 
accurate atomic calculations for analysis, and these cal- 
culations will require contributions from quantum elec- 
trodynamics effects. 

There are accurate calculations and measurements 
available for thallium QilUllHII], and measurements have 
also been considered for the Ba + ion [li| [23[ and are in 
progress for the Ra + ion [2J|. 



II. QED CORRECTIONS 

The quantum electrodynamics corrections considered 
in this work arise from vacuum polarization and elec- 
tron self-energy. The inclusion of vacuum polarization 
is numerically relatively simple, achieved via inclusion of 
the Uehling potential. The self-energy contribution is in- 
cluded via the radiative potential method developed in 
Ref. 14j. Note that the radiative correction to the elec- 
tric dipole transition operator cLei (vertex correction) is 
very small and may be neglected. The change of the 



electric dipole matrix elements come from the QED cor- 
rections to the electron wave functions. 

The parity nonconservation amplitude of a transition 
(a — b) between states of the same parity can be expressed 
via the sum over all possible intermediate opposite parity 
states n, 



EpNC = 



(b\d.Ei\n){n\hw\a) (b\hw\n)(n\dEi\a) 



E a — E n 



■ + ■ 



Eh — E„ 



], (1) 



where Aei is the electric dipole transition operator and 
h\y is the nuclear-spin-independent weak interaction. 
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TABLE I: Values (from Ref's [lOHjJ]) for the percentage con- 
tributions of vacuum polarization including the Uehling and 
smaller Wichmann-Kroll (W-K) potentials, and self-energy- 
vertex (SE-V) to weak s-p matrix elements for various atoms. 
Uncertainty is estimated from the spread of values of different 





Uehling 


SE-V 


W-K 


Total 


Rb 


0.20 


-0.51 


0.001 


-0.31(2) 


Cs 


0.40 


-0.84 


0.003 


-0.43(3) 


Ba+ 


0.41 


-0.86 


0.003 


-0.45(3) 


Tl 


0.93 


-1.44 


0.015 


-0.50(5) 


Fr 


1.13 


-1.75 


0.02 


-0.60(7) 


Ra+ 


1.17 


-1.80 


0.02 


-0.61(7) 



We consider QED corrections from the three sources 
separately - corrections to the weak amplitudes, the 
dipole amplitudes and the energy denominators. 

In the original works, QED corrections were only cal- 
culated for the weak matrix elements (see e.g. |10l - [l3| ). 
This was a reasonable approximation numerically for ce- 
sium due to a chance cancellation between QED contri- 
butions coming from corrections to the energy levels and 
dipole matrix elements, which were calculated by Flam- 
baum and Ginges This cancellation is not guaran- 
teed and it was demonstrated in Ref. [14[ that corrections 
from all three sources are equally important and must be 
included. Full determinations of QED corrections to the 
entire PNC amplitude have only been considered for ce- 
sium ninm. 

The radiative potential method, developed by Flam- 
baum and Ginges [l4|, defines an approximate potential 
L such that the radiative correction to the energies coin- 
cides with its average value, 

SE n = (n\L\n). 

This potential takes into account the local vacuum po- 
larization operator (including the lowest-order in Za 
Uehling potential as well as the higher order Wichmann- 
Kroll potential) and the non-local strong Coulomb field 
electron self-energy operator. 

In this work, we use the existing calculations of QED 
corrections to the weak matrix elements, which are pre- 
sented in Table [fl These values have been taken from the 
works Ref's Iol-[l2l] (see also [HI])- Midpoints and uncer- 
tainties have been chosen to agree with different previ- 
ous determinations. Note also that these calculations are 
valid only for s-pi/2 weak transitions. 

We then use the radiative potential method outlined 
above, with the exception of the small Wichmann-Kroll 
term, to calculate the QED corrections to the energy lev- 
els and electric dipole matrix elements. We then calculate 
the dominating terms in equation ([T]) with and without 
the radiative corrections to determine the total percent- 
age correction to the PNC amplitudes for several transi- 
tions in Rb, Cs, Ba+, Tl, Fr and Ra+. 



III. CALCULATIONS 

We use the sum-over-states approach to calculate 
the PNC amplitudes. Relativistic Hartree-Fock (RHF) 
method is used to construct single-electron orbitals and 
random phase approximation (RPA) is used to include 
the effect of core polarization by external fields. 

Core- valence correlations are included by means of the 
correlation potential method (CPM) 25]. The second- 
order correlation potential E^ 2 ) is calculated using many- 
body perturbation theory and then used to construct 
the so-called Brueckner orbitals (BO) for the external 
electron. BO are found by solving the Hartree-Fock-like 
equations with an extra operator E: 



(2) 



where Hq is the relativistic Hartree-Fock Hamiltonian 
and index n denotes valence states. The BO ip^ ^ and 
energy E n include correlations. 

We then use a simple scaling procedure to estimate 
the contribution of higher-order correlations by includ- 
ing a factor in front of the correlation potential, \YS 2 \ 
which is chosen to reproduce energy levels of the lowest 
lying valence states. A separate A is used for each of the 
initial/final states, and another is used for each set of 
intermediate states (e.g. np l / 2 , np 3 / 2 ). For the second- 
order correlation potential, values for the fitting param- 
eter typically take values A ~ 0.8 — 0.9. This fitting gen- 
erally increases the accuracy of the wave-functions and 
therefore the matrix elements. 

The PNC amplitude is given by the expression similar 
to ([1]), in which states a, b, n are single-electron RHF 
states, and operators cIei and hw are modified to in- 
clude the effect of core polarization: d,Ei — s> d,Ei + $Vei, 
hw —> hw + SVw- Here SVf is the correction to the 
self-consistent core potential due to the effect of external 
field / (/ is either electric field of laser light <1ei or weak 
electron- nucleus interaction hw)- The corrections to the 
core potential are found by solving self-consistently the 
RPA equations for the core states 



(3) 



Here Hq is the RHF Hamiltonian, index c numerates core 
states, / is the operator of external field (weak or electric 
dipole). 

QED corrections are included by adding the radiative 
potential to the RHF Hamiltonian Ho. This is done on 
the stages of calculating single-electron RHF and Brueck- 
ner orbitals and solving the RPA equations for the elec- 
tric dipole field. We remind the reader that we don't 
include QED corrections for the weak interaction. The 
more reliable results for weak matrix elements are found 
in different approaches considered before in Refs. (Iol - [l3| . 

B-spline technique [26[ is used to construct sets of 
single-electron orbitals used in the calculate E^ 2 ) and for 
summation in ([T]). The basis states used to calculate 
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FIG. 1: Sample contribution to the Qp\/2 — 6p 3 /2 PNC transi- 
tion in thallium arising from the double core polarization by 
weak (dashed-dot line) and electric dipole (dots line) fields. 
The 6s state is treated as a core state. 



£( 2 ) are the linear combinations of B-splines which are 
the eigenstates of the RHF Hamiltonian Hq (RHF or- 
bitals). The basis states used in the summation ([T]) are 
the eigenstates of the Ho + Hamiltonian (Brueckner 
orbitals). 

The approach described above does not take into ac- 
count the effect of core polarization due to simultaneous 
action of the weak and electromagnetic interactions. An 
example of the contribution of this type is presented on 
Fig. [TJ for the case of thallium. The contribution of the 
double core polarization to all PNC amplitudes consid- 
ered in present paper except thallium does not exceed one 
per cent. Therefore, it can be neglected in calculating the 
relative effect of QED corrections. The contribution of 
the double core polarization to the 6p 1 / 2 — 6p3/2 PNC 
amplitude of thallium as about 30% (Fig. Q]) if thallium 
is treated as a single-valence-electron atom, leaving the 
6s electrons in the core. Note that the double core po- 
larization was taken into account in our old calculations 
of the PNC in thallium [H[ by means of the RPA-like 
approach with two operators of external field. Similar 
calculations in present work would not be practical since 
we want to separate the effect of QED corrections on 
weak and electromagnetic matrix elements and it would 
be hard to do so in the RPA-like approach with two op- 
erators. Therefore, we use the configuration interaction 
(CI) method instead. The calculations for Tl are simi- 
lar to our recent calculations of the EDM enhancement 
factor for thallium [27[ . The Tl atom is treated as a three- 
valence-electrons systems. The calculations are done in 
the V N ~ 3 approximation (27|. The correlations between 
6s and 6p electrons are included very accurately in the 
framework of the CI method. The core-valence correla- 
tions are also included by means of the many-body per- 
turbation theory. The QED corrections are included in 
electric dipole matrix elements while weak matrix ele- 
ments are kept unaffected. 

Our calculations, which include core-polarization and 
(fitted) second-order Brueckner-type electron correla- 
tions, agree with existing calculations and experimental 



TABLE II: Calculations of PNC amplitudes for several tran- 
sitions in the RPA and fitted second-order Brueckner (BO) 
approximations and comparison with previous calculations. 
Values given in units i(—Qyy/N) x 10 -11 a.u. . 



Transition 


This work 
(RPA + BO) 


Previous 
calculations 


85 Rb 




0.139 


0.139 [6] 




5s-4d 3 / 2 


0.449 




133 Cs 


6s-7s 


0.897 


0.8977 [5] 




6s-5d 3 / 2 


3.75 


3.76 [19] 


138 Ba+ 


6s-7s 


0.671 






6s-5d 3/2 


2.36 


2.34 [19] 


223 p r 


7s-8s 


15.4 


15.41 [20] 




7s-6d 3/2 


59.4 


59.5 [19] 


226 Ra+ 


7s-8s 


10.9 






7s-6d 3/2 


46.6 


45.9 [19] 


values to 


around 1% for 


energy levels and 


a few percent 



for El matrix elements. Importantly, we calculate the 
PNC amplitudes which also agree well with previous cal- 
culations. A summary of this comparison is presented in 
Table HH 

IV. RESULTS AND DISCUSSION 

Percentage QED corrections to individual lowest en- 
ergy levels are presented in Table IIII) and corrections 
Rei to dipole matrix elements are presented in Table 
HVl The factor Rei is defined 

(b\d E i\a) = (6|4i|a) (l + -Rei), (4) 

where the subscript indicates the zeroth order matrix 
element, without radiative corrections, and a is the fine 
structure constant. 

Tables Mil and IIVI show QED corrections to the lowest 
states only. However, we include QED corrections to 
all states used in the summation ([1]). The results are 
presented in Table fVl together with the correction arising 
from the QED correction to the weal matrix elements. 
The latter are taken from previous works [Tol - [l^ | (see 
also Table HI). 

While there is typically still some cancellation between 
the contributions from the energy levels and dipole ma- 
trix elements in the transitions studied, it is not always 
as complete as it is with cesium, making these results 
significant. 

Despite the fact that there are some individual correc- 
tions in the atoms that are quite large, there is significant 
cancellation between contributions from the weak matrix 
elements and the combined contributions from the ener- 
gies and dipoles in most of the transitions presented. This 
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TABLE III: Percentage QED corrections to ionization ener- 
gies of lowest states for several atoms. 



Atom 


Correction (%) 


Rb 


5s 
-0.040 


6s 
-0.023 


5Pl/2 

0.003 


6P1/2 

0.002 


7Pl/2 

0.002 


0.003 


5P3/2 

0.001 


Ps 


6s 
-0.069 


7s 
-0.040 


6P1/2 

0.006 


7pi/2 

0.004 


8pi/ 2 
0.003 


9P1/2 

0.002 


5d 3 / 2 
0.031 


Rn + 
Da 


6s 
-u.uoo 


7s 

-U.UoO 


6P1/2 

u.uuo 


7Pl/2 
U.UU4: 


8pi/ 2 
u.uuo 


9P1/2 

u.uuz 


5^3/2 
U.UZo 


TV 


6^1/2 

-0.01 


6P3/2 

-0.02 


6s 
-0.14 


7s 
-0.07 








Fr 


7s 
-0.142 


8s 
-0.076 


7pi/2 

0.002 


8pi/ 2 
0.001 


6^3/2 

0.067 


7P3/2 

-0.003 


8P3/2 

-0.002 


Ra+ 


7s 
-0.109 


8s 
-0.066 


7pi/2 

0.001 


8P1/2 

0.000 


6^3/2 

0.059 


7P3/2 

-0.005 


8P3/2 

-0.004 



causes the total QED contributions in most atoms to be 
highly suppressed. 



The uncertainty of the total QED corrections to the 
PNC amplitudes come from uncertainties of all three 
sources (Hw, E n and El) added in quadrature. The 
uncertainty of the first term (Hw) has been estimated 
m 

.10 12]-' The estimation of uncertainties for two other 
terms (E n and El) comes from the spread of values of the 
QED correction found in different approximations (RPA, 
Brueckner, etc.). 

TABLE V: QED corrections (as percentages) to PNC am- 
plitudes for several atoms. Corrections due to weak matrix 
elements (Hw) are taken from Table [I] Corrections due to 
change of energy denominators (E n ) and electric dipole tran- 
sition amplitudes (El) are the result of present work. 

Transition H w E n El Total 



a V N 3 approximation 



TABLE IV: QED corrections Rei to the dipole matrix ele- 
ments of several atoms. Rei is defined in equation Q. 



Atom 


Transition 


Rb 


5s-5pi /2 
0.193 


5s-6pi /2 
-1.624 


6s-5pi/ 2 
-0.254 


6s-6pi /2 
0.180 


5P3/2-5S 

0.196 


Cs 


6s-6pi /2 
0.328 


6s-7pi /2 
-3.678 


7s-6pi /2 
-0.452 


7s-7pi/ 2 
0.301 


6P3/2-6S 
0.346 


Ba+ 


6s-6pi/ 2 
0.257 


6s-7pi /2 
11.75 


7s-6pi/ 2 
-0.544 


7s-7pi/ 2 
0.240 


6P3/2-6S 

0.283 


TP 


6pi /2 -6s 
0.619 


6pi /2 -7s 
-0.804 


6P3/2-6S 
0.717 


6P3/2-7S 

-0.527 




Fr 


7s-7p 1/2 
0.647 


7s-8pi /2 
-5.588 


8s-7p 1/2 
-0.731 


8s-8pi /2 
0.561 


7P3/2-7S 

0.768 


Ra+ 


7s-7p 1/2 
0.483 


7s-8pi /2 
26.39 


8s-7p 1/2 
-0.857 


8s-8pi /2 
0.441 


7p 3 /2-7s 
0.622 



a V N 3 approximation 



Rb 


5s-6s 


-0.31 


-0.25 


0.31 


-0.25(4) 




5s-4d 3/2 


-0.31 


0.12 


0.001 


-0.19(5) 


Cs 


6s-7s 


-0.43 


-0.42 


0.52 


-0.33(4) 




6s-5d 3 / 2 


-0.43 


0.20 


-0.003 


-0.23(7) 


Ba+ 


6s-7s 


-0.45 


-0.54 


0.68 


-0.31(4) 




6s-5d 3/2 


-0.45 


0.29 


-0.05 


-0.22(8) 


Tl 


6pi/2-6p 3 /2 


-0.50 


0.07 


0.06 


-0.37(8) 


Fr 


7s-8s 


-0.60 


-0.83 


1.02 


-0.41(8) 




7s-6rf 3/2 


-0.60 


0.42 


-0.02 


-0.2(1) 


Ra+ 


7s-8s 


-0.61 


-0.97 


1.20 


-0.38(9) 




7s-6d 3/2 


-0.61 


0.47 


-0.09 


-0.2(1) 
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